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The field equations of Mannheim's theory of conformal 
gravity with dynamic mass generation are solved numerically 
in the interior and exterior regions of a model spherically sym- 
metric sun with matched boundary conditions at the surface. 
The model consists of a generic fermion field inside the sun, 
and a scalar Higgs field in both the interior and exterior re- 
gions. From the conformal geodesic equations it is shown how 
an asymptotic gradient in the Higgs field causes an anoma- 
lous radial acceleration in qualitative agreement with that ob- 
served on the Pioneer 10/11, Galileo, and Ulysses spacecraft. 
At the same time the standard solar system tests of general 
relativity are preserved within the limits of observation. 

PACS: 04.25.Dm,95.35.-Hd 



I. INTRODUCTION 

It has become increasingly probable in recent times 
that general relativity is on its way to being falsified as 
a classical theory of gravitation. General relativity fails 
to account for the motions of galactic clusters ( Q , (|] ) 
and rotations of individual galaxies ( 0, see also 
1^ and references therein) without the ad hoc assump- 
tion of varying amounts and distributions of cold dark 
matter . At the same time the existence of cold dark 
matter in the universe is becoming more and more in 
doubt. As determined by the Boomerang Experiment, 
the small amplitude of the second peak relative to the 
first in the power spectrum of the cosmic microwave back- 
ground does not support the cold dark matter hypothesis 
1^. Even within the solar system general relativity has 
met with a potential failure. Radio metric data from 
the Pioneer 10/11, Galileo, and Ulysses spacecraft indi- 
cate a constant (with respect to radius from the sun), 
radial acceleration directed toward the sun of magni- 
tude ~ 8.5 X 10~^° m s~^, in addition to the expected 
Newtonian-Einsteinian acceleration that varies as the in- 
verse radius squared 0. Furthermore, Viking Lander 
ranging data limit any such anomalous radial accelera- 
tions on the planets Earth and Mars to less than one hun- 
dredth of that detected on the spacecraft, a violation of 
the equivalence principle. On a more fundamental plane, 
general relativity is not compatible with dynamical mass 
generation in that the latter leads to a manifestly trace- 
less energy-momentum tensor as the source for Einstein's 
field equation, while the Einstein tensor, the left side of 



the field equation, is not traceless. 

Collectively these problems pose a challenge to gen- 
eral relativity as a classical theory of gravity, and it may 
be appropriate to seriously consider alternatives. Two 
such alternatives have appeared in the literature in re- 
cent years: (1) Milgrom's modified Newtonian dynamics 
(MONO) and (2) the Weyl fourth-order, confor- 

mally invariant gravitational theory which has been re- 
examined recently by Mannheim and Kazanas ||l0[| and 
extended by Mannheim to include dynamic mass gen- 
eration Both of these theories have been used to 
account for galactic rotation curves without dark matter 
[|l2|,^. However, only (2) is compatible with dynam- 
ical mass generation, in fact requires it, and therefore we 
believe (2) is the stronger candidate as an alternative to 
general relativity. 

In this paper we investigate conformal gravity with 
dynamical mass generation in the solar system to see 
whether the theory can account for the anomalous, con- 
stant radial acceleration on the spacecraft while also pass- 
ing the standard solar system tests, namely light bending 
at the limb of the sun and perihelion precession. As the 
field equations of conformal gravity are nonlinear, fourth- 
order differential equations, exact analytic solutions are 
not possible, even in a simple solar system model. We 
first consider approximate solutions and then confirm the 
general validity of the approximate results with complete 
numerical solutions of the interior and exterior problems 
with matched boundary conditions at the limb of the sun. 
We find that a completely conformal theory, including a 
scalar Higgs field for dynamical mass generation, does, 
in fact, predict an anomalous, nearly-constant radial ac- 
celeration toward the sun on a test particle such as a 
spacecraft, while at the same time satisfying the stan- 
dard solar system tests. Furthermore, we give a qual- 
itative argument why planets should behave differently 
than spacecraft and not exhibit the same anomalous ra- 
dial acceleration. 



II. CONFORMAL GRAVITY WITH DYNAMIC 
MASS GENERATION 

Conformal gravity is based on the Weyl geometry ac- 
tion. 
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2a / ct^xyf^ -Rp^i?^" - -R^ + (surface term) , (1) 



where C^^\p is the Weyl conformal tensor, g = det {g^v), 
Rfj.iy\p is the Riemann tensor, i?^^ = R'^ is the Ricci 
tensor, R = R^ ^ is the Ricci scalar, and a is a dimension- 
less constant. The Weyl action is invariant under a con- 
formal transformation, g^^, [x) — > [Q. [x)] g^^ {x), which 
is an arbitrary, continuous deformation of the spacetime 
manifold. We use a metric signature (— , -|-, and the 
sign conventions of Weinberg . 

For dynamic mass generation, following Mannheim 
pl| , we choose a conformally invariant matter action in a 
model consisting of a self-interacting scalar Higgs field S 
and a generic fermion field if) with a Yukawa interaction 
with the scalar field as in the Standard Model of par- 
ticle physics. The requirement of conformal invariance 
restricts such a matter action to the form 

/™ = - / d^x^ [7^ {x) (d^ + {x)) + ihs] ^ 



+ ld,Sd^S-l(^^y' + XS^^ , (2) 

where (x) is a spin connection, /i is a dimensionless 
Yukawa coupling constant, and the quadratic term in S 
with i?/6 < creates a situation of spontaneous symme- 
try breaking and consequent dynamical mass generation. 

Setting the variation of the total action, I = 1^ + Im, 
with respect to the metric equal to zero yields the Weyl 
field equations. 



(3) 



where the source, 



T^, = SI„JSg>''' = 1^7^ (x) [d, + (x)] ^ 

2 1 ^1 1 ^ 



-g-S^ (Rfj.iy - -^9tivR ) - gt^u^S^ , 



(4) 



is the energy-momentum tensor, which in this theory in- 
cludes the scalar Higgs field, and the tensor on the left 
side, 

1 2 



2 11 

+ -^RRfiu — '^Rti^Rxu + -;j^9fj.vR\pR^^ ~ q9iiuR^ (5) 

was first obtained by Bach [isl] . 

Setting the variation of the action with respect to the 
two fields to zero yields the two field equations of motion. 



(x) [d^+r^ {x)]-hS}ij = Q, 



(6) 
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(7) 



Equation (|^) is the Dirac equation in curved spacetime 
with a fermion mass due to spontaneous symmetry break- 
ing given by 



m = hS. 



(8) 



Equation is the Klein-Gordon equation in curved 
spacetime with a self interaction, —AXS^, and a fermion 
source, —hipij). The scalar curvature, or more precisely 
-~R/6, plays the role of the mass squared of the Higgs 
boson. In contrast to the Standard Model of particle 
physics, where the field point of the Higgs vacuum is 
a universal constant over all spacetime, in the present 
context the Higgs field is dynamically connected to the 
metric of spacetime through Eqs. (||), and (0), and 
is not constant in general. In the present problem the 
Higgs field varies radially from the sun, and it is the ra- 
dial gradient of the Higgs field that is responsible for the 
anomalous radial acceleration. 

The traceless property of the energy-momentum tensor 
with dynamic mass generation is most easily established 
from Eq. (|^) starting in a gauge where the scalar Higgs 
field takes on a constant, nonzero value Sq. In this case, 
using the field equations of motion, Eqs. (||) and (Q), 
we see that T^^ = . Since T^^ is a conformal ten- 
sor, under conformal transformation, — g^T^^ 
[n^ (x)] g^^O-^ (x) = r^u^ ^ general 

the scalar Higgs field transform as S (x) ^ S (x) /il{x). 
Therefore, in a general gauge in which the Higgs field 
is not constant, S (x) = SQ/fl{x), and the trace of 
the energy-momentum tensor is still zero since T'^^ = 

{x) T^f, = 0. 

The Bach tensor, the left side of the field equations, 
Eq. (||), is also traceless, as can be seen from Eq. (|) 
which gives 



2 a 2 

^R' ;A+3i?-% 



2R^^:a: 



l-L,X 



- 2R^'^R 



= 0, 



where we have used the covariant divergence of the 
contracted Bianchi identity, {R''^ - gf'^R/2) = 0, 
to establish the cancellation of the third and fourth 
terms. Thus dynamical mass generation is compati- 
ble with conformal gravity but not with general relativ- 
ity since the Einstein tensor is not traceless in general: 



lg^^^R = -R^O 



In general relativity the predictive power of the theory 
lies in the postulate that a test particle of mass m follows 
a world line on the spacetime manifold that minimizes the 
action [1161 



dr . 



(9) 
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where dr = ^ — g f^^dx^^ dx^ is the hne element. With 
the mass m of the test particle constant, the path is a 
geodesic. However, the action of Eq. (||) is not suitable 
for a conformal theory because it is not invariant un- 
der conformal transformation. In fact, conformal gravity 
has been criticized by Perlick and Xu as being non- 
predictive except for null geodesies because the action of 
Eq. is arbitrarily variable under conformal transfor- 
mation. What these authors have failed to recognize is 
that a conformal theory must be completely conformal. 
The action of Eq. must be replaced by a confor- 
mally invariant action, and the way to make this change 
is apparent from Eqs. (||) and (^. The mass of the 
test particle in a conformal theory must no longer be a 
fixed attribute of the particle, but instead it is a dynamic 
quantity generated through spontaneous broken symme- 
try and given by Eq. (||). The mass of the test particle 
now depends on position in spacetime and must be moved 
inside the integral. 



elm [x] dr = —he I S {x) dr 



(10) 



The action of Eq. ( p^ is invariant since under con- 
formal transformation dr ^ f2 [x) dr and S {x) — > 
[17(a;)]"'S'(a;). Thus we see that a conformal theory 
of gravity is not only compatible with dynamical mass 
generation but is non-predictive without it. 

Conformal geodesic equations are obtained from the 
conformal action given by Eq. ( p!o| ) by requiring it to be 
stationary for arbitrary small variations in path between 
two fixed points in spacetime: 



5S_ 
he 



S {x) dr 



Pi 



[SSdT + SS (dr)] = , 



Pi 



Pi 



dx"^ 



dx^ dx'^ dg^u 
dr dr dx'^ 



dxf'diSx'') d{5x'')dx'' 

+9t^<y— 7Z ^ 9ay 



dr dr 



dr dr 



dT = 0, 



and, integrating by parts with Sx"' (Pi) = 5x'^ {P2) — 0, 
we have 



1 d 



Pi 
Pi 

S 



dS 1 ^dxi" dx" dgf,^ 
dx'^ 2 dr dr dx'^ 



dx" 

d7 



Sx^dT = . 



For arbitrary variations Sx"^ and multiplying by g'"^ / S 
and summing on cr, we obtain the conformal geodesic 
equation, 



d'^xP (b^dx^ 1 dS 



dr dr 



(11) 



In their original paper on conformal gravity Mannheim 
and Kazanas obtained an analytic solution of the 
homogeneous equations. 



W^. = , 



(12) 



for a static, spherically-symmetric metric corresponding 
to the line element. 



ds^ (r) e^dt^ + -r^dr^ + (d0^ + sin^ Odcf)^) 

b(r) ' 



It is given by 



5(0 = 1-^^^^-3/^,+ 



— nr , 



(13) 



(14) 



where /3, 7, and k are integration constants. 

As Mannheim and Kazanas have shown, the above 
simple form of the metric suffices in a conformal the- 
ory. Starting with the general line element in terms of a 
radial coordinate p. 



ds^ = -B{p) e^dt^ + A (p) dr^ 



(15) 



one makes a change of radial variable p — p{r) and 
rewrites the line element as 



ds' 



p{r) 



1 2 



r 

„2 /jq2 



-b (r) e'^dt'^ + a (r) dr^ 

-sin^ edcj)^)] , (16) 



where 6 (r) = r'^ B {p (r)) / [p (r)]^ and a (r) — 
r'^A{p{r)){dp/dry / [p{r)]. Then, if one chooses the 
fmiction p (r) to satisfy 

If dr , , 

(17) 



p{r) 



r^^A{p{r))B{p{r)) ' 



the line element becomes 



ds' 



'pir)' 


2 r 


r 





-b (r) e^dt^ + 



b{r) 



dr' 



W (d6'2 +sin2 6ld02)] , 



(18) 



which is related to the simple line element of Eq. ( |l3| ) by 
a conformal transformation with 17 (r) = r/p (r). 

Mannheim and Kazanas [0 interpret Eq. (^4|) as an 
exterior vacuum solution with a Schwarzschild-like met- 
ric in a background de Sitter spacetime, where the in- 
tegration constant 7 measures the departure from the 
Schwarzschild metric. They have argued that with 7 <C 1 
the metric of Eq. ( p^ ) enjoys all of the successes of gen- 
eral relativity in the solar system, while on larger in- 
terstellar scales in a galaxy the linear term ^r can pro- 
vide the increasing gravitational potential with radius 
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required to supplement the decreasing Newtonian poten- 
tial in order to account for the plateau characteristics 
of galactic rotation curves They have also argued 

that the integration constant k is so small that the de 
Sitter term is negligible except on cosmological distance 
scales. Perlick and Xu [0 correctly have insisted that 
the values of the integration constants, (3, 7, and k, can- 
not be chosen at will, but must be determined by match- 
ing boundary conditions with an interior solution at the 
inner boundary with the source and at infinity. 

But there is an even more serious criticism of the vac- 
uum solution of Mannheim and Kazanas, and that is it is 
not a vacuum solution at all. The vacuum in conformal 
gravity is not T^^ = because even outside the source 
the Higgs field must be nonzero. Otherwise a massive test 
particle would not be possible. Even in a gauge where the 
Higgs field is constant, outside the source the final term 
in Eq. (^) contributes to a nonzero energy-momentum 
tensor. The next to last term also contributes for all met- 
rics in which the Ricci tensor is nonzero, including the 
metric of Mannheim and Kazanas given by Eq. . In 
a general gauge all the terms in Eq. except the first 
contribute outside the source. We have argued above 
that conformal gravity not only is compatible with dy- 
namical mass generation through spontaneous symmetry 
breaking of a Higgs vacuum, but also requires it in order 
to be predictive. It is therefore not logically consistent to 
ignore the Higgs contribution to the energy-momentum 
tensor outside the source, and the vacuum solution of 
Mannheim and Kazanas is not physically relevant by it- 
self. All we can say is that mathematically it is an exact 
complimentary solution of the homogeneous field equa- 
tions that may be a good approximation to an exterior 
particular solution. 



tensor will be small outside the sun if the radial gradient 
of the Higgs field S and the self coupling constant A are 
both small with respect to unity. We will see later from 
our numerical calculations that both of these assump- 
tions are valid in the asymptotic region. Note that for the 
Schwarzschild metric the next to last term in T^^ is iden- 
tically zero since the scalar curvature R = i?^^ = 0. But 
realistically the metric must be near to Schwarzschild, 
but not exactly Schwarzschild because the mass squared 
of the Higgs particle is i?/6. 

Assuming spherical symmetry such that S — S (r) and 
for a line element of the form of Eq. (nsh with 



b{r)^l 



2m 



(19) 



where m = GM/c^ is the geometric radius of the sun, 
the conformal geodesic equation, Eq. (|ll]), evaluated for 
p = 0, 1, 2, 3 is given respectively by (c = 1) 



IK 



= 0, 



(20) 



2 



dtV 
d^J 



1 db 
26 



br 



or sm I — 

dr 



1 dS 
S'd^ 



0, (21) 



d^ 



2 dr de 
r dr dr 



— sm cos ( 



dr 



1 dS de dr 
S dr dr dr 



= 0, 



(22) 



III. APPROXIMATE ANALYTIC 
CONSIDERATIONS 

Before presenting our matched interior-exterior numer- 
ical solutions we begin our investigation by asking the 
question, what are the characteristics that would be re- 
quired of the Higgs field in order to produce an anoma- 
lous, constant radial acceleration at large radii in the 
solar system such as that observed on Pioneer 10/11 etc. 
spacecraft, while at the same time permitting the the- 
ory to pass the standard tests in the inner solar system. 
This question can be answered by the conformal geodesic 
equation. We wish initially to answer it analytically, so 
we have to make a few simplifying assumptions. First, 
since the theory must pass the standard solar system 
tests, the metric must be very close to the Schwarzschild 
metric, and we will assume initially that it is exactly 
Schwarzschild. One way to assess this assumption is that 



the Schwarzschild metric satisfies Ru 



0, and therefore 



it also satisfies W^i, =0. So the Schwarzschild met- 
ric is a first approximation to the conformal metric with 
~ 0. As we see from Eq. the energy momentum 



dT^ 



2 cot 6*- 



de , 



dr dr 



2 dr d4> 
r dr dr 



1 dS d(j) dr 
S dr dr dr 



0. (23) 



For spherical symmetry, we lose no generality by re- 
stricting motion to the = 7r/2 plane, in which case Eq. 
( p^ ) is trivially satisfied. Equations ( po| ) and ( p3| ) may 
be integrated once {I and k are constants of integration) , 



dt 
d7 



I 

bS 



2 k 

and the result substituted into Eq. ( ^l| ) giving 

2' 



(24) 



(25) 



}_dS_ 



1 db 
26 dr 



52 



dr 
d7 

P6 



(26) 
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From the line element with d9 = and l/bS substituted 
for dr/dr according to Eq. (|2j), we write 



1 = 



bS^ 



(27) 



Solving Eq. 
gives 



7^ for [dr/dr] and substituting into Eq. 



}_dS_ . 2 
^3 dr ^ 



and finally for radial motion (fc — 0), 



dr^ 



1 db P dS 
2d^ ^ 53 d7 



0. 



(28) 



(29) 



As the radiometric measurements of the accelerations 
of the spacecraft in the solar system are based on co- 
ordinate rather than proper time, we convert Eq. ( |29| ) 
to coordinate time using dt/dr = l/bS and restore c by 
letting t ^ ct, 



di^ 



1 db 
b dr 



1 dS 
S'd^ 



b'^c^S^ db b^c^ dS 



2P dr S dr 



= 0. 



(30) 



It is the final term in Eq. ( |3(]| ) that has the potential 
of accounting for a constant radial acceleration for large 
radius. In this limit & = 1 — 2m/ r — > 1, so we set 



, din S* 
dr 



a(r) 



(31) 



where a (r) is an anomalous acceleration. The next to 
final term in Eq. (^0|) accounts for the Newtonian accel- 
eration. 



b^^S^ db 1 



2P dr 



2m 



GM _^ S'2 GM 



(32) 



since for r = 40 AU, 2m/r = 2.0 x 10~^ < 1. Turning 
our attention to the second term in Eq. ( |30| ) , we sec from 
Eqs. (|l|) and ^ that 



1 db 

b dr 



ldS_ 
S dr 



a{r) 



^ 2 m 



2GM , A /Idr 



which is negligible for nonrelativistic speeds. 

Based upon the above stated assumptions and approx- 
imations, the conformal geodesic equation, for nonrela- 
tivistic radial motion with a Schwarzschild metric, re- 
duces asymptotically to 



df^ 



52 GM 



(34) 



The requirement, based upon the radiometric data from 
the spacecraft, that a (r) uq be asymptotically con- 
stant determines the asymptotic form that the scalar 
Higgs field must take through Eq. (|3l|): 

d In 5 ao 



dr 



— InS* = + ln5o , 5 = 5*0 exp ^-|r^ 



(35) 



For ao = 8.5 x 10~^° m s"^ and r = 40 AU, aor/c^ = 
5.6 X 10~^^ <C 1 and the required asymptotic form is well 
approximated by 



S" 5*0 ( 1 + -^r 



Then the Newtonian term in Eq. (| 
totically 



(36) 

becomes asymp- 



GM 
'p' — 



2 GM 



GM 



(37) 



if we set the two integration constants equal, I = Sq, and 
neglect a^r /(? with respect to 1. Note from Eq. ( p4| ) that 
I sets the scale of proper time, and we are free to measure 
proper time in any units we wish. Finally, if the Higgs 
field takes on the asymptotic form of Eq. (|35|), then the 
conformal geodesic equation becomes in the asymptotic 
region 



di2" 



GM 



(38) 



While Eq. ( pq ) gives the required asymptotic form of 
the Higgs field in order for the conformal theory to ac- 
count for the anomalous radial acceleration, it remains 
to be determined whether the conformal field equations 
have such an asymptotic solution. The answer to this 
question requires the full interior and exterior solutions 
with matched boundary conditions at the limb of the 
sun. But before turning our attention to a full solution 
of the field equations, it is instructive to consider an ap- 
proximate analytic solution of Eq. (Q), the equation of 
motion of the Higgs field, in the exterior region (-0 = 0), 
assuming the metric is exactly Schwarzschild (i? — 0) and 
neglecting the self-interaction term (A = 0). In this case 
and with spherical symmetry Eq. (|^) becomes 



'^'^ dr^ 



(f^^dS_ 
dr ' 



2m 



2mr 



r . 



(33) ^^^^ differential equation has a solution. 



(39) 



(40) 



Although there are two constants of integration in the 
solution of the second-order differential equation, physi- 
cally only their ratio is important, as can be seen from 



S^Ci\n{l- — \+C: 
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the conformal geodesic equation, Eq. (|1 1|) , in which for 
the present static case with spherical symmetry the Higgs 
dependent factor in Eq. is 



ldS_ 
S dr 



2m 



1 - 



2m 



In 1 - 



2m 



(41) 



Thus we can choose either Ci or C2 arbitrarily and it is 
convenient to evaluate C2 by setting 5 = 5*0 = 1 m^^ at 
r = i? at the limb of the sun. Then we have 



Co 



1 - Ci In 1 



2m 



(42) 



Ideally we should evaluate Ci by matching the slope, 

2mCi 



dS_ 

dr 



2mr ' 



(43) 



with the slope of an interior solution at the limb of the 
sun. But since at the present stage of the exposition we 
do not have an interior solution, we instead evaluate Ci 
by the observational requirement of an anomalous radial 
acceleration of ao = —8.5 x 10^^" ms^^ at a radius of 



i?4o = 40 AU, giving 



Ci — 5*0 



2 m 



2mc 
aoR^ 



■40 



2m Y 



(44) 



Evaluating Eqs. 



and (E3), we obtain Ci 



1.20 X 10-'' m-i and C2 = 1 + 5.14 x 10" 



1 



With these constants, the Higgs field S is found to have a 
value that is nearly constant for r > 0.1 AU, with a slope 
that is positive and very small with respect to unity (see 
Fig. |l|, noting that only the first half AU is shown for 
clarity, which establishes the asymptotic trend). This ap- 
proximate solution of the scalar field equation of motion 
shows that the features the scalar field can be expected 
to have are intuitively reasonable. Although the slope 
is not constant as would be required by observation, the 
inclusion of a nearly, but not exactly, Schwarzschild met- 
ric and a small but nonzero self coupling is enough to 
reduce the variation with radius significantly, as we shall 
see. The approximations used (zero Higgs mass and self 
coupling), although quite sweeping, are not so extreme 
as to take the analytic solution of the reduced equation 
far from the solutions of the full equation. 



IV. SOLAR SYSTEM TESTS 

Two further questions need to be considered, and these 
are the effect of the scalar field on the deflection of light 



and perihelion precession, the standard solar system tests 
of any gravitational theory. These questions can be an- 
swered by the orbit equation (u = 1/r and h is the an- 
gular momentum per unit mass), 



d2 



u db 1 



S^dh 



bS 



dS 
du 



(45) 



which is derived from the conformal geodesic equation, 
Eq. (pi]), for the case of the conformal line element given 
by Eq. (p^). Regarding the deflection of light, we note 
that the angular momentum per unit mass of a photon 
is infinite, and therefore the final term in Eq. ( |45| ) in- 
volving the scalar field is zero. This result is compatible 
with the fact that the photon, being massless, does not 
couple directly with the Higgs field. Thus there is no 
direct effect on the defiection of light. Indirectly the pre- 
diction of the conformal theory will differ from that of 
general relativity by the degree that b (r) differs from the 
Schwarzschild metric. We will see in the numerical so- 
lution that this difference is minimal. On the question 
of perihelion precession, while the final term in Eq. ( |45| ) 
is nonzero we note that it does not contain the inverse 
radius, u. In the standard iterative solution of Eq. ( |45| ) 
assuming small eccentricity, the zeroth order approxima- 
tion, u = uo (1 + e cos (/)), is substituted into the right- 
hand side of the equation, resulting in a secular term, 
(3m^//i'*) ecj) sin (j), that grows with the angle traversed. 
Again, the presence of the Higgs field does not contribute 
directly to the secular term and consequently to perihe- 
lion precession. Only through the deviation of b (r) from 
Schwarzschild will the prediction of the conformal theory 
differ from general relativity. 



V. NUMERICAL SOLUTION 

Finally we present a numerical solution of the Weyl 
field equations and the scalar field equation of motion 
that has been carried out over both the interior and ex- 
terior domains, with boundary condition matching at the 
limb of the source, a prototypical sun which is presumed 
static and spherically symmetric. The exterior domain 
of the numerical solution has been extended in fine reso- 
lution to a radius of 60 AU, and in coarser resolution to 
further radius. The simple solar model used to perform 
these solutions was a polytrope of order three, and the 
equations were solved for the metric coefficient b and the 
scalar field S, with the metric coefficient a removed by 
coordinate and conformal transformations. The fermion 
field was handled by following the averaging procedure 
used in ||lll| and solving in terms of generic fermion num- 
ber density and pressure. The main methods of solution 
were Runge-Kutta integrations and Newton-Raphson re- 
laxations. 

Boundary conditions at the limb of the source were 
provided by requiring compatibility with observation. 
For the metric field, the Schwarzschild solution is known 
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to be very accurate in the inner solar system for light de- 
flection and precession observations, so the metric field 
was required to match well with the Schwarzschild solu- 
tion at the limb of the sun. The scalar field has it's value 
normalised to one at this point. The limb gradient of the 
scalar field can be limited in magnitude by compatibil- 
ity with existing observations, or a more specific value 
can be produced by reference to the anomalous motion 
of the Pioneer spacecraft. If the observed acceleration 
of these spacecraft is accurate (and a transponder signal 
should be the most sensitive test we have of these ef- 
fects to date), the scalar field gradient corresponding to 
the observed acceleration at 40 AU would be S' ~ 10^^'' 
m^^. The exterior numerical solution gives a correspond- 
ing limb value of S' ~ 10""'^^ m~^. 

With these boundary conditions in hand, a consistent 
interior /exterior numerical solution was produced (Fig- 
ures |2[ P). The interior metric field merges well into an 
exterior solution of the Schwarzschild form. The scalar 
field increases across the interior domain before smoothly 
assuming a very slight gradient as required in the exte- 
rior region. Newton-Raphson iterations on the exterior 
solution near the source show that the metric field pro- 
duced is in agreement with Schwarzschild to one part in 
10^^, even in the presence of a scalar field. From our 
arguments given above, we would conclude that theory 
passes the standard solar system tests to this accuracy. 

The far-field (30 - 60 AU) behaviour of the fields is 
also satisfactory. The metric field continues to agree with 
Schwarzschild to good order, up to the accuracy limits 
of the numerical methods. The scalar field assumes a 
very slight and asymptotically decreasing gradient which 
varies very slowly across the exterior domain. The re- 
sults are very similar to the restricted analytic solution 
described above, but with the inclusion of the metric field 
in the mass term resulting in a slower variation of the 
scalar field gradient. It is expected that with improving 
accuracy of the far-field methods that the scalar field will 
vary more slowly again. 

The values that the constants of the theory assumed 
with these boundary conditions were a ~ 3.5 x 10^^ and 
A < lO^'*^. The value of a is set through the boundary 
conditions, and that of A by the stability of the solution at 
large radius, larger values producing runaway behaviour 
in the fields at large radius due to the XS^ term in the 
equations. 

VI. PLANETARY MOTION 

The form we have found for the scalar field shows a 
radially and rapidly decreasing gradient in the vicinity of 
the source. This feature could have relevance as to why 
anomalous accelerations are not observed for planetary 
bodies, at least to the accuracy of ranging experiments 
to date ( 0; see also comments in [|l8|). While it is not 
yet known how the scalar field sourced by more than one 



body would be configured, we expect that a secondary 
body, such as a planet, would superpose a similar char- 
acteristic well upon the background field produced by the 
primary. Such a superposition may well result in an av- 
erage of the scalar field gradient over the volume of the 
planet that is less than the single-source field gradient. 
Thus the constituent particles of the planet could, on av- 
erage, experience less of an anomalous radial acceleration 
towards the sun than they would if the rest of the planet 
were not present. In this respect a planet is unlike a 
spacecraft, which would not appreciably affect the scalar 
field in its own vicinity. 

VII. CONCLUSION 

In this paper we have investigated the theory of con- 
formal gravity with dynamical mass generation, as origi- 
nally presented by Mannheim The theory has been 
presented here in a logically consistent format that clar- 
ifies some areas in the construction of the theory that 
have been confused in the past. In particular it has 
been emphasized here that since conformal gravity not 
only is compatible with dynamical mass generation, but 
also requires it to be predictive, the scalar Higgs field 
is an integral part of the theory and cannot be ignored. 
From this framework, analytic and numerical results have 
been extracted. The import of these results is that it ap- 
pears that conformal gravity with dynamical mass gener- 
ation can reproduce the gravitational effects within the 
solar system, including the possible Pioneer spacecraft 
accelerations, while still satisfying the standard solar sys- 
tem tests. This successful outcome was made possible 
through the inclusion of the scalar field required for dy- 
namical mass generation. Solutions of the full system of 
equations have not been obtained before. 
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FIG. 2. Scalar Field - Boundary Matching 
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FIG. 3. Metric Coefficient - Boundary Matching 



